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^ ■ Abstract 

We investigate the nonperturbative behaviour of the quark propagator in 
£S) | axial gauge using a truncation of the Dyson- Schwinger equations which re- 

spects the Ward-Takahashi identity and multiplicative renormalisability. We 
demonstrate that above a critical coupling a c , which depends on the form of 
the gluon propagator, one can obtain massive solutions for both explicit and 
dynamically generated quark masses. The stability of these is discussed in 
the context of the Cornwall-Jackiw-Tomboulis effective action. 
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I. INTRODUCTION 



ON 

^3 
D., 

The Dyson-Schwinger equations (DSEs), in effect the functional Euler-Lagrange equa- 
tions of quantum field theory, provide a natural framework for investigating nonperturbative 
Green functions. Whilst this infinite tower of equations encodes all information relating 
Green functions of different order, their use was restricted until quite recently because of 
the often brutal approximations which were necessary to render them soluble. A great deal 
of progress has been made subsequently, much of which is reviewed in Ref . |1|] , which com- 
prehensively covers QED, QCD and applications to hadronic physics. We will be concerned 
with the DSE for the inverse quark propagator, shown diagrammaticaliy in Fig.|I|, which 
relates the fully dressed quark and gluon propagators to the complete quark-gluon vertex. 
Early studies replaced the gluon propagator and quark-gluon vertex by their perturbative 
expressions, violating the the Ward-Takahashi identity (WTI) connecting the quark prop- 
agator to the vertex and destroying the universality of the QCD coupling constant. More 
recently, Bashir and Pennington, extending previous work by Curtis and Pennington, have 
shown how imposing the constraint that the critical coupling marking the onset of chiral 
symmetry breaking be gauge covariant places strong restrictions on the transverse compo- 
nent of the quark-gluon vertex @,[|. The longitudinal part can of course be determined 
from the WTI ||. Whilst studies of DSEs in such a scheme appear to be at an early stage 
due to the complex nature of the expressions involved, the weaker condition of requiring 
that multiplicative renormalisability be respected, which also puts strong constraints on 
the transverse vertex, has been used widely and with considerable success; particularly in 
massive quenched strongly-coupled QED 4 ||. 
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Lattice simulations hold out the possibility of a first principles determination of non- 
perturbative propagators ||. However, their use has been hampered because with current 
technology they are unreliable below ~ 1 GeV, where the most interesting and novel be- 
haviour is expected to lie. Additionally, gauge fixing on the lattice reduces the number of 
gauge configurations available for analysing whilst at the same time introducing ambiguities 
due to the presence of Gribov copies. The latter could be avoided by implementing an axial 
gauge fixing procedure on the lattice, rather than the more usual Landau gauge, but unfor- 
tunately there are then further complications associated with imposing periodic boundary 
conditions. It would seem therefore that at the present time DSEs are the most reliable tool 
for studying the infrared behaviour of propagators in the continuum limit. 

In what follows we consider the quark propagator S(p), given in terms of functions F{p) 
and M{p) by 

and with the form for the quark-gluon vertex specified by Curtis and Pennington ||. In 
axial gauge the propagator may have additional structure of the the type --y.n^.pG + H), n 
being the four- vector appearing in the axial gauge fixing condition n^A^ = 0. The addition 
of such terms, although potentially important, results in a large increase in the complexity 
of the problem, and represents a direction for future research. We choose also to work with 
p.n = as a further condition on the specification of the gauge, again for practical reasons to 
make the equations more tractable. The general gauge dependence of the quark propagator 
is extremely interesting, but lies beyond the scope of present DSE studies except in the case 
of QED 4 in covariant gauges, where details are now beginning to emerge j7],|| 



II. THE QUARK DSE 

The quark DSE for S(p) at spacelike momenta, — p 2 > 0, may be written formally as 

l=rpS(p) ~^-J D flu (k-p)S(k)T»(k 1 p)S(p)d% (2) 

where a is the strong coupling, C/ is the quark Casimir, D^k — p) is the full gluon prop- 
agator and Y u (k,p) is the quark-gluon vertex. The ansatz which we shall use for the latter 
function is 

V(k,p) = Tl(k,p)+TKk,p), (3) 

with 

rvjfe v ) = "L ( J— + 1 ^ ■ i(k+ P y( r k + rP ) ( i i ^ 



2 \F{k) F(p) J 2 k 2 -p 2 \F(k) F(p) 

( M(k) _ M(p) \ (k + pY 
\F{k) ~ F{p) J k 2 -p 2 



(M(k) M(p) yk+ P Y (1) 



and 
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in which D(k,p) is a regular function of and p which behaves as k 2 for k 2 ^> p 2 . The 
transverse contribution IV cannot be specified uniquely. However, its form may be strongly 
constrained by demanding that the quark propagator respects the requirements of mul- 
tiplicative renormalisability, as described by Curtis and Pennington. We have described 
previously the insensitivity of solutions of the quark DSE to the precise form of D(k,p) 
and will adopt here 

k 2 + p 2 

which guarantees that multiplicative renormalisability is respected to all orders in leading- 
and next-to-leading logarithms. This is the form suggested by Curtis and Pennington in 
their important study of the structure of the transverse vertex and is the most widely used 
in recent DSE work. Upon substituting expressions (||) for the vertex into the quark DSE, 
it is clear that in the absence of any mass terms, i.e. with M(p) = 0, there is a single linear 
integral equation for the wavefunction F(p) which looks generically like 

j = F(p) _ p j* ^pW^(MFW z {{k _ p?) d% (7) 

where the shorthand 47T/9 = aCf has been used and the kernels Ai(k,p) and A2(k,p) 
are functions to be specified. Z(k 2 ) represents the nonperturbative aspect of the gluon 
propagator, as will be described in more detail in the next section. The form of equation 
(0) follows from noting that the integral on the right-hand side of the DSE is proportional 
to F(p)F(k), whilst the vertex contains terms which are inversely proportional to F(p) or 
F(k) separately. A is an ultraviolet regulator in the guise of a sharp momentum cutoff. 

In earlier work we elected to evaluate the angular integrals analytically by means of a 
shift in variables to remove the angular dependence from the gluon propagator and then by 
assuming that F(p) was sufficiently flat that the approximation F ((k — p) 2 ) « F(k 2 + p 2 ) 
could be made |J. This idea had been used previously in a study of the gluon propagator 



in axial gauge ||10|| , but its validity is highly questionable and we shall show that it leads to 
spurious errors in the solution for F(p). 

We have indicated elsewhere how unique, finite solutions to (|7|) exist provided that the 
function <fi(p) defined by 

m = 1 ~Pj ^%^Z{k 2 ) d'k (8) 

is nonzero for all p 2 ||. Since A 1 (k,p) and therefore 4>(p) is independent of F(p), this 
implies the existence of a value of (3 for which can be made zero and which results in 
the breakdown of solutions for F(p) with M(p) = 0. If the integral in <p(p) is positive, 
this (3 defines a critical value of the strong coupling a which is interpreted as signalling 
chiral symmetry breaking. When the mass function M(p) is reintroduced, <ft(p) is changed 
to become 
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which has an integrand suppressed at the origin relative to <p(p) provided Mik 2 — > 0) ^ 0. 
This lifts the condition <p(p) = and means that solutions for F can be recovered beyond 
the critical coupling exhibited above. M(p) may be introduced by hand as an explicit quark 
mass, or it may be generated by the dynamics of the quark-gluon interactions. 

The F{p) appearing in equation ([?]) is a bare quantity which has an inherent dependence 
on the ultraviolet regulator A 2 and should be written strictly as F(p 2 ;A 2 ). Its subtrac- 
tive renormalisation follows from introducing a renormalisation constant F^fj, 2 , A 2 ) so that 
F(p 2 ; fi 2 ) = F^lfi 2 , A 2 )F(p 2 ; A 2 ) is the renormalised wavefunction which depends on p 2 and 
a renormalisation point fi 2 at which we impose the boundary condition F(/j?; fi 2 ) = 1, cor- 
responding to the free propagator. The numerical solution of the resulting integral equation 
for F(p) is amenable to a number of methods. The one used here is to map the integration 
over k onto a finite logarithmic mesh, which entails both infrared and ultraviolet cutoffs 
whose effect will need to be determined. Following this the integral is replaced by a discrete 
sum so that 

N 

/ K(p, k)f(k) dk^J2 K (Pi> k j)f(kj) Uj (10) 

where the weights Uj are chosen according to a Gaussian quadrature rule and A" is the 
number of grid points used. 

At each value of (p, k) it is necessary to have also a numerical evaluation of the angular 
integrals since the gluon contains the dependence Z ((k — p) 2 ) = Zik 2 + p 2 — 2pkcos6). 
With these steps, the solution of (|7|) becomes equivalent to determining the elements f(pi) 
which satisfy A" simultaneous linear equations of the general type 

TV 

f( Pi ) + £ K( Ph k s )f{k s ) u>j = g( Pi ). (11) 

3=1 

The accuracy of the discretised approximation can be found retrospectively by iterating the 
equation for F(p) with increasing A^ to ensure that the error is small. 



III. CHIRAL SOLUTIONS 

The kernels Ai(p, k) and A 2 (p, k) do not depend on either F(p) or Z{k 2 ) and are given in 
the appendix. An expansion of Ai and A 2 for k 2 ^> p 2 reveals that the integrand involving 
Ai decreases as a power {k 2 )~ s Z {k 2 ) with 5 = 3/2 and the integrals involving it are finite. 
On the other hand A 2 behaves as (A; 2 ) for k 2 ^> p 2 and the finiteness of its contribution 
depends on the asymptotic behaviour of F(k) and Z(k 2 ). We must specify the precise form 
of Z(k 2 ) to be used to determine the critical coupling and numerical solutions for F(p). An 
additional point worth noting is that in a scheme where MR is not enforced, the divergence 
in the equations is present in the Ai rather than the A 2 terms. The coefficient of the 
divergence is a function of p 2 and cannot be removed by a straightforward subtraction. 
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In one of the first serious studies of the gluon DSE, Baker, Ball and Zachariasen (BBZ) 
examined the consequences of a gluon propagator which had the same tensor structure as the 
perturbative one combined with the scalar function Z(k 2 ) representing the nonperturbative 



dynamics so that [|TT | 



2 . Z(k ) / k^n u + kyTt^ 2 k^k v 



c ^ ) = _^^_=_^ +n ^. (12) 

More generally in axial gauge the propagator could contain the additional contribution 

^ 2 )=P 2 (^)(^-^), (13) 

where p2{k 2 ) is a further function to be determined. The value of taking ( |I2D alone is that 
the terms in the gluon DSE involving the full four-gluon vertex are projected out, resulting in 
a considerable simplification. Of course, the main attraction of axial gauge is the absence of 
the ghost fields, which removes one of the diagrams involved in the gluon equation and means 
that the Ward- Slavnov- Taylor identity for the quark-gluon vertex can be replaced by the 
Ward-Takahashi identity, enabling us to use the expression for the vertex in ([5]), which was 
originally proposed for QED, without further approximation. The solutions found by BBZ, 
which incorporated only the longitudinal contribution to the quark-gluon vertex, exhibited a 
Z(k 2 ) which diverged strongly as (/c 2 ) -1 in the infrared and which had the asymptotic limit 
(In fc 2 )~ 11//16 , which disagrees with renormalisation group improved one-loop perturbation 
theory. However, in attempts to describe phenomenology by the exchange of one or more 
soft gluons, infrared divergences arise unless the gluons are softer than (k 2 )' 1 as k 2 — > 0. On 
the other hand, the axial gauge propagator cannot be nonsingular at k 2 = 0. The remaining 
possibility is for a cut solution which behaves as (k 2 )'^ with < ( < 1. The nonlinearity 
of the gluon DSE suggests that there may exist alternative solutions even within the same 
truncation schemes. 

A very successful phenomenology of diffractive scattering, elastic scattering, the gluon 
structure function and total cross-sections results from the idea, first suggested by Land- 



shoff and Nachtmann [12], of a process- independent softening of the gluon propagator in the 
infrared domain jn|. Models of diffractive scattering, such as two-gluon exchange, suffer 
from infrared divergences if the propagator is not softer than 1/A; 2 as k 2 —>■ 0. In particu- 
lar, it is not possible to obtain the experimentally observed slope of the hadronic scattering 
cross-sections in the t-channel with an extremely infrared-divergent gluon. Whilst a strongly 
infrared- singular propagator might well be relevant to the question of the nature of confine- 
ment, a softer behaviour is far more phenomenologically efficacious. 

With this is mind, we wish to examine the consequences of a gluon propagator which is 
less singular than the one suggested by BBZ and which may therefore be phenomenologi- 
cally relevant. Several alternatives are considered. As indicated above, it is believed that 
the gluonic wavefunction cannot be singular at the origin and the perturbative (quenched) 
propagator with Z a (k 2 ) = 1 thus provides one limiting case. At the other extreme, one 
could have a propagator with no singularity at the origin, i.e. with Z(k 2 ) = k 2 as k 2 — > 0. 
Of course, the gluon cannot behave like this in the ultraviolet regime and it is necessary to 
introduce a scale /j, 2 which marks the transition from the nonperturbative to the perturba- 
tive region. Perturbatively, Z(k 2 ) runs as [lnfc 2 ] -1 asymptotically. In order to reproduce 



5 



this behaviour and to have Z(k 2 = fi 2 ) equal for all choices of the propagator, we take as 
our second ansatz 



^"{Snfe-l + x 2 )]- 1 , x 2 >l (M) 

with x 2 fi 2 = k 2 . 

Z a and Zfj represent two extreme alternatives and an intermediate between them is re- 
quired. Ideally this should share the correct UV behaviour of Z& and lie between Z a and 
for k 2 < fi 2 . One possibility, which does have these properties, has phenomeno logically desir- 
able infrared behaviour and which can be applied successfully to reproduce the experimental 
data for p — p scattering is that proposed by Cudell and Ross [M : 



x 2 

Zc ^ X ^ = 0.88(x 2 ) - 22 - 0.59(x 2 ) - 86 + 0.95a; 2 ln(2.1a; 2 + 4.1) ' ^ 

This was derived originally from a study of the nonlinearity of the gluon DSE, which admits 
the existence of more than one solution. That particular analysis is flawed because of a sign 
error [ 18 1 , but the phenomenologically attractive properties of equation (pl|) plus its position 
between Z a and Zb in the infrared domain make it an ideal ansatz for present puurposes. 

The broad shape of the curves obtained for F(p) is not affected by Z(k 2 ). They are 
illustrated in Fig. |2| with the choice Z(k 2 ) = Z c (k 2 ), for a = 0.2,0.6, 1.0, 1.4. These show 
a plateau for p 2 < fi 2 , whose magnitude increases with a, passing through the fixed point 
at F(fi 2 ) = 1 and thereafter decreasing logarithmically as (lnp 2 ) -5 with 5 ~ 0.2. The 
magnitude of F(p) at large p 2 is suppressed as a is increased. In fact the value of a when 
F(p) crosses the axis is coincident with the loss of solutions and corresponds to the critical 
coupling. In this case the critical value a c was found to be a c ~ 1.45. For the alternatives 
for Z(k 2 ), Z a and Zb with the infrared limits Z a — > 1 and — > k 2 , we found a c = 0.9 and 
a c = 1.9 respectively. Interestingly the former of these is close to the a c ~ 0.93 which has 
been found in Landau gauge investigations of the quark propagator with a quenched gluon. 
When previously we used the approximation F((p — k) 2 ) « F(p 2 + k 2 ) to derive the integral 
equation for F, we found oscillations in the region p 2 ~ fi 2 , whose magnitude grew with the 
coupling a until they dominated the solutions at a ~ a c . They are completely absent from 
the present set of solutions and are thus clearly an artifact of the approximation. This is 
not too great a surprise, for whilst the derivative of F(p) is indeed small for p 2 -C fi 2 and 
p 2 ^ fi 2 , F{p) changes rapidly near p 2 = fi 2 . 

The value of a c is dominated by the infrared behaviour of the gluon and reduces with 
an increasingly singular gluon, which is an interesting effect in its own right. It suggests 
that the more rapidly the gluonic wavefunction increases at large distances, the lower the 
coupling at which chiral symmetry is dynamically broken. A determination as to whether 
the behaviour of the gluon is the major factor contributing to chiral symmetry breaking 
would require a more detailed analysis in which the quark and gluon DSEs were solved 
simultaneously. Conversely, it is quite plausible that an even softer gluon would not support 
chiral symmetry breaking at all. The test of whether the breakdown of chiral symmetry is 
the correct interpretation of a c will follow from examining the more general case of a massive 
propagator. 
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IV. THE MASSIVE DSE 



Apart from the obvious complication of adding mass terms to the quark DSE, a feature 
which has been neglected frequently is the correct renormalisation prescription in the pres- 
ence of a chiral symmetry breaking mass. This is nontrivial and of great importance to a 
fully self-consistent analysis of the solutions for F(p) and M(p). It is worth presenting the 
details here and we follow the notation used in Ref. ||. In what follows, tildes denote renor- 
malised quantities and primes denote regularised ones which depend on the UV regulator 
A. 

The starting point is the renormalised inverse quark propagator S^Qo) which can be 
written in terms of the renormalised current mass m(/z) and renormalised self-energy function 
E(p; //) as 

S _1 (p) = >y.p - m{p) - E(p;/i) = [y.p - m (A)] Z 2 (ji, A) - E'(p; //, A) (16) 

where all the dependencies have been shown explicitly, Z%(ji, A) is the fermion renormalisa- 
tion constant and mo (A) is the bare mass, which is zero in the massless limit of QCD. 
The self-energy can be decomposed into Dirac and scalar components as E(p;/x, A) = 
7.pEi(p; A) + E 2 (p;/i, A) with a similar expression for E'. At the renormalisation point 
the boundary condition 

S _1 (p) =j.p-m(jj.) (17) 

is imposed. Upon inserting the decomposed self-energy into and using (pD, a compar- 
ison of the scalar terms and those involving 7 matrices reveals the relations between the 
renormalisation constant Z2, Ei, E2 and the bare and renormalised masses to be 

Z 2 (ii,A) = l + Ei(/z;^,A), 
Z 2 {/i, A)mo(A) = m(/i) - E' 2 {/i; /1, A). (18) 

One can see from the second of these that in general for finite A, m(/i) = does not neces- 
sarily imply that Z 2 (fi, A)m (A) = - this only happens when m(/z) = E^/x; //, A). As the 
cutoff is removed however it is necessary that limA^oo Z 2 mo = to ensure conservation of 
the axial current jg whose divergence is d^j^ = d^ipj^j^ip)- 111 practice we are forced to 
adopt a finite cutoff in numerical simulations, but one can check the dependence of Z 2 iuq 
on A to verify that it does tend to zero in the limit. 

Using (|16D together with the relations between bare and renormalised quantities in ([18]) 
it is straightforward to show that the relation between the renormalised and regularised 
self-energies is 

E(p;/x) = E , (p;/i,A)-E , (//;//,A). (19) 

Writing the quark DSE in terms of renormalised quantities and recalling that in axial gauge 
Z\ = Z 2 we can see that the regularised self-energy is given by 

S'(p;^A) = zZ 2 (/i,A)0i J K d A kYmf°(k,p)D pa (k-p). (20) 
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Taking the trace of both sides, or the trace after multiplying throughout by ^.p yields a pair 
of coupled nonlinear equations for the two components of We wish to find the propagator 
in the form specified by (|1|), in terms of which the renormalisation condition for M(p) is 

M(p;fi) M(p;/i,A) M(/i;/x,A) 

— = 171(11) + — -. rv — rr. (21) 

F(p;/x) ^ F>(p;fi,A) F'(//;/x,A) V ; 



A. Explicit Quark Mass 



Explicit chiral symmetry breaking (ECSB) can be induced by introducing a non-zero 
renormalised quark current mass m(/x). Performing the traces in the DSE produces an 
equation for F(p) of the form 



A[^,k) 



dk = F(p) - (3F(p) J 



Ai(p,A0 



k 2 + M 2 (k) ^' ^ ^' J k 2 + M 2 (k) 

A' 2 {p,k) - A' 2 (/i, k) 



dk 



F(k) dk, 



k 2 + M 2 (k) 

which remains linear in F, and a nonlinear equation for M of the form 

m( \ (\v(\<R7( a^/ ^ f ^ u Xi(p,k)M{p) + X 2{p,k)M{k) 
M(p) = m{fi)F{p) + (3Z 2 {fi,A)F{p) J dk — - 

R7 l ^ [ ^ u Xi(v,k)m( y fi) + X2(v,k)M( y k) 
- [3Z 2 ^, A)F(p) J d k P + M2(fc) 



(22) 



(23) 



The kernels A' 1; A' 2 ,Xi and \2 have been relegated once more to the appendix and we will 
refer to the denominator k 2 + M 2 (k) as Di(k, M) or simply as D\. A complete determination 
of the quark propagator now entails the simultaneous solution of the two integral equations 
( ^2|) and (23). The one for F is simpler as it is still linear and can be solved, given a 
particular M, using the method outlined for the massless case. Equation (p3|) for M, being 



nonlinear, is more difficult and the approach taken was to introduce a generalisation in the 
form of a parameter A into the denominator Di which now becomes Di = k 2 + XM 2 (k). The 
advantage of this is that A = defines a linearised equation which is easy to solve and we 
can then increment A from to 1, solving iteratively at each stage, to obtain the required 
solutions. The equations for F and M are solved alternately so that a given F^ is used 
to obtain an approximation M^' which is put back into the equation for F to find F^ t+1 ^ 
etc. until convergence of both functions is achieved. The criterion used was that the error 
function defined by 
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j 3 



3 I 3 



(24) 



satisfied £j < 10" 
for M. 



for all points pj at which F was determined, with an equivalent condition 



S 



As an illustration of the functions obtained, the equations were solved for p 2 = 1 with 
m(p) = 0.2 using the Cudell-Ross form for the gluon propagator. The wavefunction F and 
mass function M are shown in Figs.|3] and |] respectively for a = 0.4, 0.8, 1.2, 1.6, 2.0. The 
cutoffs used were A 2 R = 10~ 6 and Ay V = 10 6 (the infrared cutoff was not expected to be 
important because all the integrands are strongly suppressed near k 2 = and indeed there 
was no significant variation in any of the quantities calculated for A 2 R = 10~ 6 , 10 -8 , 10~ 10 
or 10~ 12 ). 

The first point of note is that as advertised, incorporating M(p) does allow us to recover 
solutions for F(p) for a > a c ~ 1.45. The magnitude of M(jp) as p 2 — > increases as a is 
increased, whilst the tail decreases with increasing a and becomes close to zero, although it 
appears never to cross the axis. Second, for a given a < a c , F{p) is both infrared suppressed 
and ultraviolet enhanced relative to the corresponding massless solution. Third, the large-p 2 
variation of F(p) with respect to a is reduced relative to the M = solutions. 

In order to investigate the dependence of the solutions on the UV cutoff, the equations 
were re-solved with the same values of p 2 and m(p) for each of log 10 A 2 = 6.0, 6.6, 7.2. ..12.0 
and the values of M(p 2 -> 0),M(p 2 = 10 6 ),F(p 2 0) and F(p 2 = 10 6 ) established. These 
were used to find mean values and deviations which are shown in Fig.[5| for a ranging from 
to 2.0, where the error bars represent ±1 standard deviation. The mean value of M(10 6 ) 
is the only one to show appreciable variation, for a > 1.5. To eliminate the possibility of 
a systematic dependence on A, the correlation p between each of the four quantities shown 
in Fig. |5| and In A were evaluated, p can take values between ±1 corresponding to perfect 
positive or negative correlations. In no case did it exceed (9(10~ 3 ), which confirms that 
there is little if any A dependence in the renormalised solutions. 



B. Dynamical Mass Generation 

The ECSB solutions discussed above where m(p) ^ are important - because quarks do 
after all have a mass - but the symmetry breaking has been put in "by hand" . To explore the 
transition between massless and massive phases, we must solve the equations for m(p) = 
(the chiral limit) to see if for a > a c a mass is generated dynamically. That is to say, can 
quark-gluon and self-interactions generate a mass from the initially massless theory? For 
m(p) =0 the equation for M becomes 

M W = P ZM A)F(p) f ^^)«MW|-^^)1*^) , (25) 

where Xi an d X2 are unchanged, as is the equation for F apart from terms involving M(p) = 
m(p) = which vanish. As noted previously, for a finite cutoff, m{p) = does not imply 
that the product Z 2 m = 0. However, X2(a*, k) is strongly suppressed for large p 2 , leading 
to only a small residual value of Z 2 mQ. This provides the basis for the success of the many 
studies which neglect the /i-dependent part of the right-hand side of (|25|) and take p 2 = Afjy. 
They incur the undesirable side effect that M(jp) is barely down from its infrared magnitude 
even at p 2 = A 2 [15], contrary to what one would expect for a dynamically generated mass 
term which should be small in the perturbative regime. 

The new equation for M was solved using a Powell hybrid method, after discretisation 
of the integrals, which is a variation of the Newton-Raphson technique which avoids the 
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evaluation of derivatives JT(|. The technique combines the advantages of a multidimensional 
root searching method with a minimisation of sum of squares to improve global convergence. 
This is convenient because it can be tailored to exclude the trivial M = solution. Solving 
the equations for M and F iteratively is problematic since for a > a c , the existence of a 
solution for F depends strongly on M. With a method which incorporates minimisation, we 
simply introduce a "penalty" to disfavour any intermediate M which makes the equation 
for F insoluble. 

Once more the renormalisation point /i 2 = 1 was used at which M(p) must be zero to 
satisfy fl2"T|). Fig.|B| shows the curves for M obtained for a ranging from 1.4 to 2.0 with 
A 2 = 10 s . It was not possible to find solutions with M^fl for a < 1.39, which is slightly 
lower than the critical a c = 1.45. For 1.39 < a < 1.45 there are two possible sets of solutions 
for F and M - one having M = and the other M/0. In Fig.|6] we see that M crosses the 
axis at p 2 = /x 2 , attaining a minimum at p 2 ~ 2 and then declining towards zero as p 2 — > A 2 . 
The size of M in the infrared shows a familiar plateau with M(0) increasing with a. This 
is consistent with preventing a zero in the function <f>{jp) governing the existence of solutions 
for F. The depth of the minimum increases with a and |M m i n | « 0.45 |M(0)|. 

The corresponding curves for F, which are not shown, are very similar to those in Fig|| 
but exhibit an interesting feature, namely that F is quite insensitive to a in the presence of a 
dynamically generated mass. It appears that the dynamical M holds F close to its massless 
form immediately below a c . As a increases, the magnitude of M increases to compensate, 
keeping <f>{jp) slightly above zero. The oscillations seen in M have also been identified in 
strongly coupled massive quenched QED 4 , where they are much smaller in size. 

The cutoff dependence of the solutions was tested for A 2 ranging between 10 6 and 10 12 and 
the renormalised M again showed no significant change. The product Z 2 m was computed 
at each value of a and A 2 and the variation in this is shown in Fig.|7|. The slow decrease is 
striking and the reason for it lies in the nature of the integrands involved in calculating 
Recalling that Z%(ji,A) = 1 + £'(//; A*, A) we have, after some rearrangement, from (EDf) 




z ^ A) = i 1+/3 i W^liWij • (26) 

with A(p,k) = A[(p, k)F(p) + A' 2 (p,k)F(k). Expanding the integrand in the case where 
k 2 ^> ji 2 we find 



A(p,k) 2 



[1 + F(k)] , (27) 



k 2 + M 2 (k) kink 2 

where the limits Q\ = 6% = and 26 2 = 9± = irZ(k 2 ) have been used for the angular integrals 
(as they appear in the appendix), together with the leading term in Z{k 2 ) ~ [/c 2 lnA; 2 ] . 
The dominant part is thus $(&) ~ [A; In A;] -1 and the contribution of this to Z 2 is 

with A 2 ^> a 2 ^> fi 2 , giving /(<&) = In In A — In In a. So ^2^0 does go to zero with increasing 
A, but only as a double logarithm - the very slow decrease shown in Fig.[7] which also shows 
a line fitting to In In A to demonstrate the close agreement. 
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V. STABILITY OF THE SOLUTIONS 



The existence of simultaneous chiral symmetry preserving and chiral symmetry breaking 
solutions presents the question as to which is the physically important stable state. This 
problem can be addressed partially through the use of the Cornwall- Jackiw-Tomboulis (CJT) 
effective action Tcjt ||17|| . The condition f§ = locates the stationary points of the CJT 
action and is in fact equivalent to the DSE for the quark propagator. In other words, 
solutions to the fermion DSE automatically give §§ = 0. At the stationary point, the CJT 
effective action can be written in terms of the propagator S(p) as Jj]] 

TcMS] = Tr [lnS- 1 ] + ^Tr [l - S^S] , (29) 

which may be evaluated for each of the sets of solutions where both M = or M ^ are 
possible. The trace indicates integration over momenta and summation over all discrete 
indices. We can consider the difference Ar between these, i.e. 

Ar = r C jT[M = o]-r C jT[M^o] (30) 

and a negative value with Ar < is interpreted as a signal that the occurrence of chiral 
symmetry breaking solutions lowers the energy of the vacuum and is preferred physically. 
Indeed, when Ar was evaluated for explicit and dynamical quark masses, it was found that 
massive solutions are always favoured. This resolves the ambiguity for 1.39 < a < 1.45, 
where both a massless quark and one with a dynamically generated mass were possible. It 
appears that the latter is favoured. 

One could go further and find the derivatives of Tcjt to ascertain the nature of the 
stationary point. However the CJT action is unsatisfactory because it is not bounded 
below. In particular one can show analytically that in Landau gauge the solutions to the 
gap equation with dynamical masses are unstable saddle points of the action although they 
still have Ar < as compared to the chiral solutions. It is not possible therefore to address 
conclusively the matter of the stability through the CJT formalism. 



VI. CONCLUDING REMARKS 

Axial gauge solutions to the quark DSE have been found, in a truncation scheme which 
respects MR, for several possible forms for the gluon propagator in the infrared domain. 
These reproduce the familiar perturbative behaviour at high momentum and exhibit a criti- 
cal value for the strong coupling a c , whose value decreases as the gluon propagator becomes 
increasingly infrared singular, independently of its ultraviolet running. Since chiral symme- 
try breaking is believed to be connected intimately with confinement, this provides further 
support, should such be required, that the confinement mechanism is linked closely to the 
infrared sector of QCD. 

Above the critical coupling, solutions to the DSE can be recovered by the introduction 
of a quark mass, either dynamically or explicitly, as has been found in numerous covariant 
gauge studies. The correct treatment of the renormalisation of the massive case has not 
always been handled rigorously in the literature, despite the efforts of several authors, and 
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yet it is crucially important in obtaining a self-consistent treatment of the DSE. It is worth 
reiterating that in numerical studies which inevitably involve a finite infrared regulator, it 
is not enough to simply set the bare and renormalised quark masses simultaneously to zero, 
although they must vanish as the cutoff is removed in order to ensure that the axial current 
is conserved. 

Although we have discussed briefly the question of stability, as indicated the CJT for- 
malism is not in itself a definitive tool. There have been a number of attempts to address 
this problem through alternative choices for the definition of the effective action and these 
represent an important direction for future investigation. 

Finally, it is known how the ansatz used for the quark-gluon vertex can be improved 
through considerations of gauge invariance. Studies incorporating a suitable vertex should 
prove extremely interesting although they will require the solution of much more complicated 
integral equations. This is particularly so in axial gauge, and will become especially acute 
when the extra structure involving 7.71 is included in the form for the quark propagator. 
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APPENDIX A: INTEGRAL EQUATION KERNELS 

Collected together here are the expressions for the various integral equation kernels which 
have been encountered earlier. The following shorthand notation has been adopted for the 
angular integrals: 



01 


= r z( P 2 +k 2 

Jo 


— 2pk cos 8) sin 2 8 cos 8 d6, 


82 


= r z( P 2 +k 2 

Jo 


-2pk cos 8) sin 2 8 dO, 


3 


= r z( P 2 +k 2 

Jo 


— 2pk cos 8) cos 8 d8, 


8a 


= f Z{p 2 + k 2 
Jo 


-2pk cos 8) dd. 



(Al) 



For the massless case with the Curtis-Pennington vertex one has 
Ai(p, k) = -2pk 2 8 3 - A—8 X - 2—8 3 + Ak 3 8 2 + Ak 3 8 A 



+ 
+ 



D(p,k) 
1 



V 



-Ap 2 k 3 8 2 + lQpk^ - 8— #1 - 4k°8 2 



p 



A 



.5 



p 



k 6 



- - -2p 6 k% + 2p 2 t% + Apk^ + Apk% - 2—8 3 + 2k b 8 2 
k z — p z I p 

A 2 (p, k) = -2pk 2 8 3 - A—8 l - 2— # 3 + Ak 3 8 2 + Ak 3 8 4 

V V 
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+ 
+ 



1 



D(p, k) 
1 

k 2 — p 2 



k 6 

Ap 2 k 3 9 2 - 16pk% + 8—^1 + Ak 5 9 2 

p 



k 6 



2 P 3 k 2 e 3 - 2 P 2 k 3 e 2 - A P k% - A P k% + 2—e 3 - 2k 5 e 2 ) 

p 



In the massive case, the kernels appearing in the equation for F are given by 

a; = -2pk 2 e 3 - A—e l - 2—e 3 + Ak 3 e 2 + Ak 3 e A 



p 



p 



+ 
+ 
+ 



i 



k 5 



D(p,k) [ 
1 

k 2 — p 2 
M 2 (k) 
k 2 — p 2 



-4 P 2 k 3 e 2 + iQ P k% - 8—^i - Ak 5 e 2 

p 



k 6 



-2p 3 k 2 6 3 + 2p 2 k 3 6 2 + Apk% + Apk% - 2—6 3 + 2k 5 6 2 

p 

- Ap 2 k9 A + Apk 2 9 3 - A—0 3 + 8k 3 9 2 + Ak 3 9 A 
p 

.k\ k 4 



A' 2 = -2pk% - A—6 1 - 2—6 3 + Ak s 6 2 + Ak 3 6 A 



p 



P 



+ 
+ 
+ 



1 



k 6 



D(p,k) 
1 

k 2 — p 2 
M(p)M(k) 
k 2 — p 2 



4 P 2 k 3 e 2 - iQ P k% + 8— + 4k 5 e 2 

p 



A 



l6 



2 P 3 k% - 2 P 2 k 3 e 2 - A P k% - A P k% + 2—e 3 - 2k 5 e 2 ) 

p 

k 4 

Ap 2 k9 A - Apk 2 9 3 + A—d 3 - 8k 3 9 2 - Ak 3 9 A 
p 



and those from the equation for M are 



Xi{P,k) 



X2(p,k) 



A P 3 k 2 e 3 + Ap z k A e 2 + Ap z k 6 e 4 + spk^ 



.2.3 



F(k) 1 
F(p) k 2 — p 2 

+ A P k% - Ak 5 e 2 - Ak b e A 

1 

k 2 — p 2 

+Ak n e 2 + Ak 5 e 4 



Ap 3 k 2 9 3 - Ap 2 k 3 9 2 - Ap z k 3 9 A - 8pfc 4 0i - Apk^9 3 



21,3/ 



+ 
+ 



1 Fjp) - F(k) 
F(p) D(p,k) 

1 Fjp) - F(k) 
F(p) k 2 — p 2 

- 2k 5 e 2 - 2k 5 e A 

i 



Ap 2 k 3 6 2 - 8pk% + 12k 5 6 2 

2 P 4 ke A + 6 P 2 k 3 e 2 + Ap 2 k 3 e A + A P k% 



+ Trfr (F(p) + F(k)) -2p 2 k6 A + Apk% + Ak 3 6 2 - 2k 3 6 A 
b (p) 
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FIGURES 




-iSp l (p) P 7" iSf(k) V(k,p) 



FIG. 1. Diagrammatic representation of the Dyson- Schwinger equation for the inverse quark 
propagator S^ 1 (p). 



1.6 




FIG. 2. Solutions to the massless quark DSE with the Curtis-Pennington vertex and Cud- 
cll-Ross gluon propagator at the couplings a = 0.2,0.6, 1.0, 1.4. The critical value in this scheme 
was found to be a c s=s 1.45. 



15 



1.8 



1.4 



0.0 



1 1 


i i i i 






X 


cc=0.4 ; 
a=0.8 - 
a=1.2 
a=1.6 ; 
a=2.0 - 


• 






— 

.... .. n\* 


- 


— ■ 






: n 2 =i 

M(n)=0.2 








»\ v 











10~ 8 10" 6 10" 4 10" 2 



10° 10 2 10 4 10 6 10 8 



FIG. 3. The wavefunction F with an explicit quark mass, m([i) = 0.2 for sub- and super-critical 
values of the coupling a. 
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FIG. 6. Dynamically generated quark mass function, i.e. with m(/x) = 0, for couplings above 
critical, showing strongly damped oscillations. 
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FIG. 7. Variation of Z 2 (fi, A)m (A) with the UV cutoff A. The dashed line is a fit to the 
predicted In In A dependence. 
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